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Evolutions of single-particle energies and Z = 64 sub-shell along the isotonic chain of = 82 are investigated in the 
density dependent relativistic Hartrce-Fock (DDRHF) theory in comparison with other commonly used mean field models 
such as Skyrme HF, Gogny HFB and density dependent relativistic Hartree model (DDRMF). The pairing is treated 
in the BCS scheme, except for Gogny HFB. It is pointed out that DDRHF reproduces well characteristic features of 
experimental Z-dependence of both spin orbital and pseudo-spin orbital splittings around the sub-shell closure Z = 64. 
Non-local exchange terms of the isoscalar a and to couplings play dominant roles in the enhancements of the spin-orbit 
splitting of proton 2d states, which is the key ingredient to give the Z = 64 sub-shell closure properly. On the other 
hand, the tt and p tensor contributions for the spin-orbit splitting cancel each other and the net effect becomes rather 
small. The enhancement of the sub-shell gaps towards Z = 64 is studied by the DDRHF, for which the local terms of 
the scalar and vector meson couplings arc found to be important. 
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During the past years the developments on the radioac- 
tive ion beam (RIB) facility have opened the new frontier 
Tor nuclear physics - EXOTIC NUCLEI, which have in- 
tensively challenged our understanding on the nature of 
nucleus from both experimental and theoretical sides. For 
such nuclear system with an extreme neutron-to-proton 
ratio, the systematics on the shell structure evolution is of 
special significance not only for the stability itself but also 
for the appropriate description of its exotic modes such 
as the halo phenomena This also provides the new 

criterion for the theoretical descriptions. 

As an example the sub-shell closure Z = 64 has been 
identified experimentally and the increasing shell gaps to- 
wards Z = 64 are also found on this sub-shell closure as 
well as relevant spin-orbit splitting 0]. In addition, the 
high-spin isomer study along the isotonic chain of N = 83 
also indicates such trend for the Z = 64 sub-shell gap, 
which increases from 2.0MeV to 2.4MeV as the proton 
number increases from Z = 61 to 65 While within the 
relativistic density functional (RDF) models, the sub-shell 
closure Z = 64 is not always reproduced properly. Even 
as the commonly used RDF model, the relativistic mean 
field (RMF) theory @, @] fails in describing the Z = 64 
sub-shell closure, which is strongly compressed by neigh- 
boring unphysical shell closure Z = 58 @, Q ■ 
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With the presence of Fock terms, the new RDF model 
- density dependent relativistic Hartree-Fock (DDRHF) 
theory has achieved many successes in describing the nu- 
clear structure properties 0, [l(| H, 11, 12, EH- Under 
the frame of DDRHF, the one-pion exchange and p ten- 
sor effects can be efficiently taken into account, which have 
brought significant improvements in describing shell struc- 
ture and corresponding evolution || dH . Especially with 
the inclusion of tensor p, the unphysical shell closures (58 
and 92), the common disease in previous RDF calculations, 
have been eliminated by DDRHF, as well as appropriate 
recovery of the Z = 64 sub-shell It is therefore an 
appropriate choice to study the enhancement concerning 
the Z = 64 sub-shell structure within the DDRHF theory. 

In DDRHF, the Hamiltonian contains the degrees of 
freedom associated with a scalar, tu vector, p vector, p 
tensor, n pseudo- vector, and photon (A) vector couplings 
[H, In this paper all the calculations are restricted 

in the spherical symmetric systems. The spherical Dirac 
Hartree-Fock equation for nucleon can be obtained as, 



dr'h(r,r')ip(r') = eil>(r), 



(1) 



where e is the single-particle energy (including the rest 
mass) and the single-particle Dirac Hamiltonian h(r, r') 
contains the kinetic energy /i kln , the direct and exchange 
part of the potential energy, respectively the local hP and 

June 3, 2009 



non-local /i E , 

h kin {r, r') = [a-p + f3M] S{r - r'), (2a) 
h B {r, r') = [E T (r)7s + E (r) + /3E s (r)] S(r - r'), (2b) 
Y G (r,r') Y F (r,r') \ 
X G (r,r>) X F {r ir ') J ' 



h E (r,r') 



(2c) 



Among the single particle Hamiltonian the local self- 
energies Eg, Eo and Et contain the contributions from 
the direct (Hartree] terms jfj, HH, 0, 17 1 and the rear- 
rangement terms [18|, whereas the non-local ones Xq, Xp, 
Yg and Yp are contributed by the exchange (Fock) terms 
19 . 14 1 . These self-energies contain effects of the occupa- 
tion probabilities in the BCS scheme. 

From Dirac Hartree- Fock equation Jl} , the contributions 
of the single particle energy E = e — M can be written as, 



Evi 



(3) 



where Skin, corresponds to the kinetic energy, E^ and 
are respectively the direct and exchange contributions 
of various meson (photon) nucleon couplings, and Er re- 
sponds for the rearrangement term. 

In this work, the calculations are performed along the 
isotonic chain of N = 82 by DDRHF with PKA1 [1], 
as compared to those by Gogny Hartree-Fock-Bogoliubov 
(HFB) with D1S 0, Skyrme Hartree-Fock (HF) with 
Sly4 [2l|, density dependent RMF (DDRMF) with PKDD 
[221 ]. In the Gogny calculations, the pairing effects in open 
shell nuclei are considered by Bogoliubov transformation 
and the same interaction (D1S) is adopted as effective pair- 
ing force. In the other calculations, the effects of pairing 
correlations are treated with BCS approximation and the 
density dependent zero range pairing force [23[ is utilized 
in pairing channel. The Hartree-Fock or Hartree Hamil- 
tonian is evaluated with the occupation probabilities in 
the BCS approximation. In DDRHF and DDRMF cal- 
culations, the pairing strength Vq is adopted as 900MeV 
as in Ref. [|| whereas in Skyrme calculations it is set to 
HOOMeV. Because of the fact that the studied isotones 
around Z = 64 lie in the stability valley, different pairing 
treatments would not bring substantial difference. 

In Fig. [1] are shown the proton (tt) single particle en- 
ergies calculated by DDRHF with PKAljl, Gogny HFB 
with D1S [li, Skyrme HF with SLy4 [ij and DDRMF 
with PKDD[23l for the isotone 146 Gd. In this work, the 
Gogny HFB single particle energies are obtained by di- 
agonalizing the HF Hamiltonian h in the HFB equation, 
which is constructed with the density matrix of the HFB 
solution. Besides the magic shell (Z = 50), DDRHF, 
Gogny and Skyrme models present distinct gap between 
the 7r2<i3/2 and ^2^5/2 states, namely the sub-shell closure 
(Z = 64). In Fig. Q]it seems that the recovery of Z = 64 
sub-shell is tightly related to the restoration of pseudo-spin 
symmetry. Above and below the sub-shell there are two 
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Figure 1: (Color online) Proton single particle energies for 146 Gd, 
calculated by DDRHF with PKA1 @], Gogny HFB with D1S [U, 
Skyrme HF with SLy4 [H and DDRMF with PKDD 



pseudo-spin doublets, respectively 7r2p {{kZs 1 /2, 7r2G? 3 / 2 }) 
and ttI/ ({7r2d 5 /2, 7iTg7/ 2 }). Consistent with the sub-shell 
recovery, well conserved pseudo-spin symmetry on the dou- 
blets 7r2j? and rrl/ are found in these three calculations. 
In contrast, an unphysically large gap between the Tr2d 5 / 2 
and 7rlg 7 /2 states is obtained by the DDRMF calculation, 
which suppresses the physical sub-shell structure Z = 64. 

To have a systematic understanding on that, Fig. [2] 
shows the pseudo-spin orbital splittings of ir2p (Fig. [2^) 
and rrl/ (Fig. [2Jd) as functions of proton number for 
N = 82 isotones. The results are calculated by DDRHF 
with PKAl, Gogny with D1S, Skyrme with Sly4 and 
DDRMF with PKDD. As the reference the data extracted 
from Ref. Q are also shown in Fig. [21 Along the isotonic 
chain, all the models present approximate degeneracy of 
the pseudo-spin doublet n2p less than lMeV energy differ- 
ence, while on the 7rl/ doublet the DDRMF calculations 
show a large violation of the pseudo-spin symmetry by the 
amount of 3MeV, which turns out to create the unphysical 
shell closure Z = 58 0, @| ■ As compared to the experimen- 
tal isotonic dependence, both DDRHF and Gogny calcula- 
tion show a kink at Z=64 for the doublet 7r2p, whereas for 
7rl/ the kink is found only in DDRHF calculations with 
PKAl. 

For the spin-orbit splitting of ir2d states, namely the 
Z = 64 sub-shell in Fig. [l] the isotonic evolutions given 
by these four models are also shown in Fig. [3K, in compari- 
son with the data 0] which represents remarkable enhance- 
ment (~1McV) towards Z = 65. Among the theoretical re- 
sults, DDRHF calculations present the strongest enhance- 
ment (~0.6MeV) as compared to Gogny and DDRMF re- 
sults (~0.2McV). Much different from experimental trend, 
the Skyrme calculations give nearly constant spin-orbit 
splitting even slightly quenching from Z = 58 to 64 
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Figure 2: (Color online) Pseudo-spin orbital splittings of pseudo- 
spin partners {t3si/2> 7r2(i3/2 } (a.) and {n2dc l /2> ""I97/2} (b) along 
the isotonic chain of TV = 82. The results are calculated by DDRHF 
with PKA1 0, Gogny HFB with D1S [H, Skyrme HF with SLy4 
[2ll | and DDRMF with PKDD [gj], in comparison with experimental 
data |3- 




Figure 3: (Color online) (a) Spin-orbit splittings Ai?2d of proton 7r2d 
states extracted from the calculations of DDRHF with PKA1 
Gogny HFB with D1S [H, Skyrme HF with SLy4 and RMF with 
PKDD [22l |. in comparison with experimental data Q; (b) Detailed 
contributions of Ai?2d from different channels given by the DDRHF 
calculations with PKA1. 



and DDRMF model gives the kink at Z = 58. Thus, 
the non-local mean fields in DDRHF and Gogny mod- 
els may play the important role to reproduce the general 
trend of the isotonic dependence of the spin-orbit split- 
tings. More specifically Fig. [SJs shows detailed contribu- 
tions of the spin-orbit splitting AE^d from different chan- 
nels given by the DDRHF calculations. In Fig. it 
is clearly seen that the enhancement of the spin orbital 
splitting is mainly determined by the non-local Fock terms 
whereas very little effect is induced by the local term, i.e., 
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Among the Fock channels, the isoscalar a and to mesons 
together with the Coulomb field, i.e., E® +uJ+A = E®+E®+ 
E^, play the dominant role in determining the enhance- 
ment of spin-orbit effects towards Z = 64. Compared to 
strong isoscalar exchange terms, the isovector ones from 
7r and p nuclcon couplings present small effects. In fact 
the proton shell evolution along the isotonic chain, consis- 
tent with proton configuration, is mainly determined by 
the proton-proton interaction, among which the isoscalar 
exchange terms contribute much stronger effects than the 
isovector ones. As compared to the T = tensor force 
3, EH, EE EE EH, the tensor effects here (T = 1) are 



fairly weak. In Fig. [5]d distinct cancelation between n and 
p (mainly tensor p) contributions is found from Z = 50 
to 58, where the valence protons will gradually occupy the 



7rlg7/2 state. Because of the tensor feature [111], the spin- 
orbital splittings are therefore enlarged by the one-pion 
exchange whereas the tensor p presents the opposite ef- 
fects. As shown in Fig. [3}d, such cancelation are always 
found along the isotonic chain of A = 82. The effects 
from the isovector channels are then strongly weakened. 
The non-local mean field, mainly the isoscalar exchange 
terms, therefore becomes specially significant to reproduce 
the enhancement of spin-orbit effect in 7r2d states of the 
N = 82 isotones. 

In addition if comparing the pseudo-spin and spin or- 
bital splittings in Fig. and Fig. [3^, one can find that the 
data represent evident consistency between the systemat- 
ics of pseudo-spin and spin orbital splittings, i.e., stronger 
spin-orbit effect accompanied by better preserved pseudo- 
spin symmetry. Among the theoretical calculations, such 
consistence between pseudo-spin degenerations and spin- 
orbit splitting is reproduced completely by DDRHF and 
partially by Gogny. 

Experimentally the Z — 64 sub-shell is determined 
by the energy difference AE^d between the ir2d 5 /2 an d 
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Figure 4: (Color online) (a) Energy difference Ai?/„j between pro- 
ton states 7rl/i 11 y 2 and ir2d 5 /2 extracted from the calculations of 
DDRHF with PKA1 [|, Gogny HFB with D1S [H, Skyrme HF 
with SLy4 and RMF with PKDD [22l |. in comparison with experi- 
mental data [3|; (b) Detailed contributions of AE^d from different 
channels given by the DDRHF calculations with PKA1. 



7rl/ixi/2 single-particle orbits 0]. In Fig. QJi are shown the 
energy differences AEhd along the N = 82 isotonic chain. 
For comparison are also shown the experimental data ex- 
tracted from Nagai et al. in Ref. [3j, which is consis- 
tent with the experimental systematics along the N = 83 
chain extracted by Odahara et al., from the analysis of 
high spin isomers [4|]. Among the theoretical calculations, 
DDRHF with PKA1 presents similar substantial enhance- 
ment (^0.5McV) as the experimental one (~1McV) to- 
ward Z=64 sub-shell. Fig. |4p presents the contributions 
of AEhd from different channels, given by the DDRHF cal- 
culations with PKA1. From Fig. Up it is found that the 
systematics of the energy difference AEhd is determined by 
the local term i?Locai while the isoscalar non-local terms 
Ef +0J+A present opposite contributions. Here the contri- 
butions from the isovector exchange terms become negli- 
gible. One may notice that the energy difference AEhd 
is given by two spin-up states, which has different physi- 
cal mechanism from the spin-orbit splitting of n2d states. 
From this point view, one can understand different effects 
of the exchange terms on the spin-orbit splittings and en- 
ergy difference AEhd as shown in Fig. [3b and Fig. 0b, 
respectively. 



In conclusion, the Z = 64 sub-shell evolution along the 
isotonic chain of N = 82 is investigated by relativistic 
and non-relativistic mean-field approaches. Mean fields 
in DDRMF and Skyrme HF are all local in the coordi- 
nate space, while those in DDRHF and Gogny HFB have 
non-local exchange contributions. The systematic trend of 
spin-orbit and pseudo-spin-orbit splitting of N = 82 iso- 
tones are also studied. The pseudo-spin symmetry is well 
reproduced by DDRHF, Gogny-HFB and Skyrme HF. In 
contrast, DDRMF produces too large splitting, leading to 
an unphysical shell closure of Z = 58. DDRHF and Gogny 
can also acount for the Z-dcpcndence of the pseudo-spin 
splitting. The empirical large increase of the spin-orbit 
splitting between n2d 3 / 2 and 7r2d 5 / 2 states towards Z=64 
sub-shell is reproduced by the theoretical results calculated 
by DDRHF, and also by Gogny but less extent. It is found 
in DDRHF that the non-local effects due to the isoscalar a 
and u) exchange terms play an essential role in determining 
the enhancements of the spin-orbit splitting of n2d states, 
which creates properly the Z — 64 sub-shell closure at the 
same time. For the spin-orbit splitting, the it and p tensor 
contributions cancel each other and the net contribution 
becomes rather small. Although there exists some dis- 
crepancy on the order of states, the enhancement of the 
experimental sub-shell gap towards Z = 64, namely the 
energy difference between the 7r2d 5 / 2 and irlhu/z states, 
is reproduced by the DDRHF model for which the local 
terms of the scalar and vector mesons are responsible. 
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